clear
clf

x = 5:.1:35;
xSample = normpdf(x, 20, 3);
plot(x, xSample)
f = 400; b = 0.1;

xTrue = 10:.1:25;
yMeas = f * b / 22 + 1;

pPrior = normpdf(xTrue, 20, 3);
pLike = normpdf(yMeas, f*b./xTrue, 0.3*ones(size(xTrue)));
pPost = pLike .* pPrior;
hold on
plot(xTrue, pPost/(0.1*sum(pPost)))
